The X-ray Spectrum of the z = 6.30 QSO SDSS J1030+0524 
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ABSTRACT 

We present a deep XMM-Newton observation of the z = 6.30 QSO SDSS J1030+0524, the 
second most distant quasar currently known. The data contain sufficient counts for spectral 
analysis, demonstrating the ability of XMM-Newton to measure X-ray spectral shapes of z ^ 6 
QSOs with integration times > lOOks. The X-ray spectrum is well fit by a power law with index 
r = 2.12 ± 0.11, an optical- X-ray spectral slope of aox — —1.80, and no absorption excess to 
the Galactic value, though our data are also consistent with a power law index in the range 
2.02 < r < 2.5 and excess absorption in the range < NH(cm^^) < 8 x 10^^. There is also a 
possible detection (~ 2cr) of FeKa emission. The X-ray properties of this QSO are, overall, similar 
to those of lower-redshift radio-quiet QSOs. This is consistent with the statement that the X-ray 
properties of radio-quiet QSOs show no evolution over Q < z < 6.3. Combined with previous 
results, this QSO appears indistinguishable in any way from lower redshift QSOs, indicating that 
QSOs comparable to those seen locally existed less than one Gyr after the Big Bang. 

Subject headings: cosmology; observations — galaxies: formation — galaxies: high-redshift — X-rays: 
galaxies — quasars: individual (SDSS J1030-I-0524) 



1. Introduction 

One of the most important goals in modern cos- 
mology is to understand the formation of galaxies 
and large-scale structures. According to the cur- 
rently prevalent theory for structure formation, 
biased hierarchical buildup within the ACDM 
framework (hereafter referred to as ACDM), the 
growth of large galaxies results from mergers be- 
tween smaller systems, with more massive galax- 
ies hosted within more massive dark matter halos, 
themselves the rarest, most biased overdensities 



in the underlying mass distribution. Observations 
of QSOs and other massive systems are thus an 
excellent way of testing structure formation mod- 
els. This is particularly apposite for QSOs at very 
high redshifts, where the added constraint of the 
youth of the cosmos allows for the most strin- 
gent tests. The discovery (Fan et al 2001, 2003) 
from the Sloan Digital Sky Survey of QSOs at 
z > 6, when the Universe was (under any cur- 
rently favored cosmology) less than one Gyr old, 
presents the opportunity to perform such tests. 



The first step is to determine the nature of the 
central engines in these QSOs, as ACDM models 
generally extrapolate the host halo mass from the 
central black hole mass (Haehnelt, Natarajan & 
Rees 1998). This is readily achieved using X-ray 
observations, which directly probe the central en- 
gines of QSOs, largely free of obscuration bias. 
Deep X-ray observations of z > 6 QSOs, to mea- 
sure their X-ray spectral shapes and luminosities, 
can therefore be used to examine the properties 
of their central black holes, and to compare to the 
X-ray properties of lower redshift QSOs to search 
for drivers behind the strong evolution of the QSO 
luminosity function (Boyle et al 2000). 

One such QSO is SDSS JlG3G-t-0524 at 10''30"27M 
+05'^24"55^0 (J2000), with a redshift of z = 6.30 
(Richards et al 2004). Deep imaging observations 
show that this QSO is not significantly lensed 
(Richards et al 2004) or beamed (Haiman & Cen 
2002), and that the rest-frame optical continuum 
shape and luminosity are probably similar to those 
of lower redshift QSOs (Fan et al 2001). The rest- 
frame UV spectrum shows an almost complete 
Gunn-Peterson trough, indicating that this QSO 
lies within the epoch of rcionization (Becker et 
al 2001). Further observations found evidence 
for supersolar metallicities (Pentericci et al 2002; 
Freudling, Corbin & Korista 2003) and indirect 
evidence for a formed host galaxy (Barkana & 
Loeb 2003). Although not detected in the sub- mm 
(Priddey et al 2003), the radio (Petric et al 2003) 
or the ROSAT all-sky survey, this QSO is de- 
tected in Chandra snapshot observations (Brandt 
et al 2002). In this letter, we present deep XMM- 
Newton observations of SDSS J1030-I-0524. We 
assume fl = 1, A = 0.7 and Hq = 70 km s~^ 
Mpc~^ 

2. Observations &: analysis 

SDSS J1030+0524 was observed continuously 
by XMM-Ncwton for 105ks in May 2003. The 
source was placed on axis, and observed by all the 
onboard instruments, although we only consider 
data from the European Photon Imaging Counter 
(EPIC) Metal Oxide Semiconductor (MOS) and 
pn cameras here. The data were reduced using the 
XMM Science Analysis System (SAS) v6.0 pack- 
age. The event lists were filtered to include only 
single and double events (patterns 0-4 for the pn 



and 0-12 for the MOS) with quality flag 0, and to 
remove time periods during which the background 
was excessively high due to proton flares. Light 
curves were extracted using the whole MOS and 
pn fields of view to check that all the flaring back- 
ground periods had been removed. The resulting 
effective exposure time was '^ 75ks. The source 
spectrum was extracted using a circular aperture 
of radius 45" centered on the source, and the back- 
ground spectrum was extracted from a contiguous 
region free from sources and with the same instru- 
mental background as our target. Additional light 
curves were extracted from the same region as the 
source spectra to search for short-term variability 
from the QSO. Spectral analysis was performed 
using the Xspec vll.3 package (Arnaud 1996), us- 
ing spectra binned to have at least 20 counts per 
bin. 

3. Results 

Before undertaking detailed analysis, we checked 
for consistency between the three EPIC detectors 
by fitting the MOS and pn spectra with single 
power-law models. The derived photon indices 
and normalizations agreed to within 0.25(t. We 
therefore proceeded to fit the three spectra simul- 
taneously. 

The X-ray spectrum of SDSS J1030-F0524 is 
presented in Figure 1. The final spectrum con- 
tained ~ 560 counts in total, where 340 were de- 
tected by the pn camera and 220 by the two MOS 
cameras. A summary of the results from spec- 
tral fitting is given in Table 1. For all fits we as- 
sumed a Galactic hydrogen column of 3.2x10^" 
cm^^ (Stark et al 1992). We first considered a 
single power law model. This gave a reasonable 
fit to the combined spectrum with a photon index 
of r = 2.12 ± 0.11. We speculate that the excess 
flux between 0.9-1.0 keV may arise from a complex 
of iron emission features (6.4-7.0 keV rest-frame), 
however the flux in this bin is only 2.1cr above 
the power law fit; we cannot therefore confirm or 
refute this without higher quality data. We also 
derived the value of aox, the slope of a nominal 
power-law between 2500A and 2KeV. Using both 
the 2500A flux and method to compute Uox de- 
scribed by Brandt et al (2002), together with our 
X-ray data, yields aox = —1.80. No short term 
variability in the X-ray emission was detected in 



any band, though given the low count rates in both 
EPIC detectors this is not unexpected. 

Previous studies (Reeves & Turner (2000), 
though see also Vignali et al (2003b)) have found 
possible evidence for increasing absorption with 
increasing redshift in QSO X-ray spectra, rising 



from 



102"cni-2 at z - 0.1 to > lO^^cm 



at 



2 > 2. It is plausible therefore that the X-ray 
emission from SDSS J1030+0524 is attenuated by 
absorption local to the QSO. A modest amount 
of excess absorption might be expected on the 
grounds that the Lya emission shows some HI ab- 
sorption (Barkana & Loeb 2003), although if the 
absorption is due to shocked material accreting 
onto the host galaxy halo it may not have sufficient 
heavy element content to produce appreciable X- 
ray absorption. We consider three power law plus 
absorption scenarios; a model where both F and 
Nh vary, a model where the power law slope is 
fixed at r = 2, and a model where the hydrogen 
column is fixed at Nh — lO^^cm^^. In the case 
where Nh is fixed the derived photon index, at 
r = 2.27 ± 0.2, lies within la of the photon index 
derived assuming zero intrinsic absorption. In the 
two cases where Nh is allowed to vary the best- 
fit values of Nh, although statistically consistent 
with zero, have errors too large to allow an ac- 
curate hydrogen column to be quoted. We can 
however explore the range of acceptable values. A 
confidence plot of F vs. hydrogen column for the 
fit in which both F and Nh are allowed to vary 
is shown in Figure 2. Based solely on this plot, 
we would quote la ranges of 2.02 < F < 2.6 and 
< A^H (cm~^) < 9 X 10^^. A hydrogen column 
above ~ 8 x lO^^cm"^ however requires F > 2.5; 
this is significantly higher than for any other ob- 
served QSO (see e.g. Reeves & Turner (2000)), 
though such steep slopes are seen in Seyferts (Wal- 
ter & Fink 1993). We therefore quote the 'best fit' 
parameters as F = 2.12±0.11, with no evidence for 
excess absorption, but also quote approximate ac- 
ceptable ranges on these two parameters that are 
consistent with our data, these are 2.02 < F < 2.5 
and < Nnicm-"^) < 8 x lO^^. 

We also explored a wider range of models for 
the X-ray emission from this QSO, though given 
the relatively low quality of our data we advocate 
caution in interpreting the following results, and 
do not present them in Table 1. A single ther- 
mal bremsstrahlung or blackbody model are both 



substantially worse fits than any of the power law 
models, and can be rejected. A power law plus 
thermal bremsstrahlung model on the other hand, 
with or without local absorption, gives a fit of 
comparable quality to the power law model but 
requires an X-ray spectral slope of F ~ 1.45; this 
is much flatter than the spectral slopes of other 
radio-quiet QSOs (Reeves & Turner 2000). While 
we cannot formally discount the possibility of sig- 
nificant thermal X-ray emission, we consider it 
unlikely. A power law x reflection model gives 
a marginally worse flt than a single power law 
model, but not sufficiently so to allow us to for- 
mally reject a model including reflection. A power 
law model with a gaussian line at an observed- 
frame energy of ~ 0.95KeV gives a somewhat 
better fit than a single power law model (with 
Xre<i/DOF== 1.4/23) but this improvement in fit 
quality is not sufficient to formally prefer this 
model over the single power law model. Further- 
more, the gaussian line properties are not well 
constrained; the only constraint is that the line 
width must be less than 600eV at 3(7 confidence. 
Overall, we conclude that the X-ray spectrum of 
SDSS J1030+0524 is best-fit by the single power 
law model, but we cannot determine the presence 
or otherwise of an FeKa line, and cannot rule out 
contributions from a thermal bremsstrahlung or 
reflected component. 

Finally, we note that the single power law 
model, whilst statistically acceptable, is not a 
'good' flt to the data, with a null hypothesis prob- 
ability of 6%. Inspection of the residuals in Figure 
1 shows that this is primarily due to several 'fea- 
tures' at observed-frame energies < IKeV, which 
are only apparent in the pn spectrum as the MOS 
spectra contain significantly fewer counts. These 
features (which include the possible FeKa line) 
have three possible origins; they could be intrin- 
sic to the QSO, they could arise from foreground 
(i.e. z < 6.3) sources, or they could be residu- 
als from (for example) soft proton flares that were 
not fully removed by the data reduction. Of these 
three reasons the flrst two are the most likely, but 
we cannot examine the reasons behind the relative 
poorness of flt further without higher quality data. 



4. Discussion 

The advent of XMM-Newton and Chandra 
have revolutionized X-ray studies of high redshift 
QSOs. Currently, detailed X-ray spectral infor- 
mation is available for a large sample of QSOs in 
the redshift range < z < 5 (Laor et al 1997; 
Reeves & Turner 2000; Reeves et al 2001; Ferrero 
& Brinkmann 2003; Vignali et al 2003b; Grupe et 
al 2004). We first compare results from these stud- 
ies to the X-ray properties of SDSS J1030+0524. 
Our best-fit spectral index, at F = 2.12 ± 0.11, 
lies near the middle of the range of values of F 
found for lower redshift radio quiet QSOs (RQQs) 
studied with ASCA (Reeves & Turner 2000), and 
is comparable to the values of F derived for z ^ 4 
RQQs from XMM-Newton observations (Ferrero 
& Brinkmann 2003; Grupe et al 2004). The 
derived rest-frame luminosity (which has a to- 
tal 1(7 uncertainty of ~ 15%) is approximately 
1.5 times and 3 times higher than those derived 
from Chandra observations by Mathur et al (2002) 
and Brandt et al (2002) respectively (who assume 
F = 2.0). We attribute this difference to the lower 
S/N of the Chandra data rather than the differ- 
ence in values of F, though variability is also a 
possibility, and one that our data do not allow 
us to test for. Our 2-lOkeV luminosity is how- 
ever comparable to lower redshift RQQs (Reeves 
& Turner 2000) and marks SDSS J1030+0524 as 
being no more than averagely luminous in the 
X-ray. Our value for the optical- X-ray spectral 
slope, aox = —1.80, is well within the observed 
range for lower-redshift RQQs, and is statisti- 
cally identical to the mean value of aox for 2: ~ 4 
RQQs (Vignali et al 2003a). Overall therefore, 
the X-ray luminosity and spectral shape of SDSS 
J1030+0524 appear to be indistinguishable from 
those of RQQs at < z < 5 (see e.g. Fig. 8 of 
Vignali et al (2003b)). Although statistics based 
on one object are obviously not trustworthy, this 
result is consistent with the statement that the 
X-ray properties of optically selected RQQs show 
no evolution up to z = 6.3. The only exceptions 
would be if (1) the unlikely event that this QSO 
has a very steep intrinsic X-ray spectrum attenu- 
ated by heavy absorption, such systems are rare 
amongst the radio-quiet QSO population, and (2) 
if the marginal excess flux in the 0.9 - l.OkeV 
energy range really is FeKa emission, as this is 
very rare in QSOs generally. Iron K emission 



from QSOs and Seyferts is thought to arise either 
from a cold reflected component from the accre- 
tion disk (for line widths ^ 50eV), or from hot 
gas in a halo (for line widths > lOOeV, in this case 
we might also expect to see significant thermal 
X-ray emission) . The detection of Iron K lines in 
SDSS J1030-h0524 would therefore be especially 
interesting, but given the quality of our data we 
do not consider this further. 

We can combine this result with previous stud- 
ies of SDSS J1030+0524 to establish whether its 
global properties differ from lower redshift RQQs. 
As described previously, this QSO does not appear 
to be significantly magnified by lensing (Richards 
et al 2004). Together with its measured X-ray lu- 
minosity and the spectral shape, this constitutes 
compelling evidence that the mass estimate of ~ 
2 X 1O^M0 for the central black hole derived from 
the 1450A magnitude (Fan et al 2001) is accu- 
rate. Considered together with previous observa- 
tions (see §1 for a review) then SDSS J1030+0524 
appears to be indistinguishable in any way from 
the lower redshift RQQ population. The very ex- 
istence of such a system only 860Myr (under our 
cosmology) after the Big Bang poses a significant 
challenge for structure formation theories, and it 
is this theme we explore in the remainder of this 
discussion. 

To explain the existence of SDSS J1030+0524 
requires the formation of a ^ IO^Mq black 
hole and probably also a ^ lO^'^Mo DM halo 
(e.g. Haehnelt, Natarajan & Rees (1998)) and 
~ 1O"M0 of stars (Magorrian et al 1998), aU 
within 860Myr. The formation of a suitably mas- 
sive halo is readily achieved within the ACDM 
framework (Mo & White 2002), and simula- 
tions predict that the mass profiles in the inner 
10/i^^kpc of the most massive halos evolve very 
little at z < 6 (Fukushige & Makino 2001; Gao et 
al 2004) . These simulations however only consider 
the hydrodynamical evolution of the dark matter 
distribution, and do not consider the astrophysical 
processes of star and black hole formation. The 
formation of the stars and central black hole must 
therefore be considered separately. There is evi- 
dence, both from observations of old ellipticals at 
z = 1.5 (Peacock et al 1998) and from the observed 
upper bound on the line-of-sight velocity disper- 
sions of stars in ellipticals (Loeb & Peebles 2003), 
that ACDM must be capable of forming the stars 



in a giant elliptical galaxy in less than one Gyr. 
This requires (for example) a 'burst' of IOOOM0 
yr~^ star formation lasting ^lOOMyr. Such high 
instantaneous star formation rates are inferred to 
exist in the 1 < z < 4 sub-mm survey sources 
(e.g. Borys et al (2003)), and in z > 4 QSOs 
(Isaak et al 2002), and lOOMyr is a reasonable up- 
per age limit for starbursts based on observations 
of local starbursts (Farrah et al 2003). To form 
the host galaxy of SDSS J1030+0524 therefore 
appears feasible, though this inference is based 
on observations of lower redshift systems, and we 
note that the upper limit on the sub-mm flux from 
SDSS J1030-h0524 (Priddey et al 2003) implies an 
upper limit on the instantaneous star formation 
rate of 3OOM0yr~^. The formation of the central 
black hole is however more difficult to explain. 
Assuming Eddington limited exponential growth, 
then a 10^ Mq black hole can grow from a IOOM0 
'seed' black hole in ~ 725Myr, however the like- 
lihood that the accretion rate is 'fine-tuned' to 
the Eddington limit for many e-foldings appears 
small, especially considering the role of feedback 
from the formation of stars in the host galaxy 
(Burkert & Silk 2001). It seems therefore that, 
unless the accretion rate exceeded the Edding- 
ton limit for some period of time, or the QSO 
luminosity currently exceeds the Eddington limit, 
accretion onto an initially stellar mass black hole 
is unlikely to produce the central black hole in 
SDSS J1030-I-0524 within the required timescale, 
and that a significant fraction of the black hole 
mass must be built via some other mechanism. 
We briefly mention two of a variety of possibilities 
(Rees 1984; Haiman & Quataert 2004) here. The 
first is that a more massive seed black hole could 
form as result of collision runaway in dense young 
star clusters (Portegies Zwart et al 2004; Gurkan, 
Freitag & Rasio 2004), which can produce 'in- 
termediate' mass black holes of several thousand 
solar masses on rapid timescales. Another pos- 
sibility is that mergers between ~ lOOM© black 
holes created in supernovae of high-mass popula- 
tion III stars contribute to the build-up of more 
massive black holes. Numerical simulations (Abel, 
Bryan & Norman 2002) suggest that one -- IOOMq 
star could form rapidly in each ~ IO^Mq ACDM 
'minihalo'. Population III stars are not expected 
to lose much mass in the final stages of stellar 
evolution and may thus be expected to form black 



holes of ^ IOOMq in less than IMyr. A fraction of 
these black holes are predicted to be driven to the 
inner regions of larger galaxies by ongoing mergers 
(Madau & Rees 2001). The space density and the 
further fate of these black holes is uncertain, but 
it is at least plausible that mergers between these 
black holes first form 'intermediate' mass black 
holes which then contribute to the build-up of the 
^ lO^M© black holes in z > 6 QSOs by a mix- 
ture of further merging and accretion (Volonteri, 
Haardt & Madau 2003). 

In summary, there is observational and theoret- 
ical evidence that the host halo, stellar mass, and 
central black hole required to make a fully formed 
QSO can form in less than a Gyr. The existence 
however of a radio quiet QSO at z — 6.3 that ap- 
pears indistinguishable from other RQQs at lower 
redshifts is still surprising, and the most press- 
ing, and as yet unanswered question is whether 
the formation of the halo, stars, and central black 
hole in an object such as SDSS J1030+0524 can 
be accomplished together in less than a Gyr in the 
ACDM framework. 
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Table 1 
Results from fitting the X-ray spectrum of SDSS J 1030+0524 
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/0.5-I2 /0.5-2 /2-IO 

(10^-^^ ergs cni"2 s"^) 
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Power law 


1.50/24 


2.12 ±0.11 
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13.70 


6.28 


6.73 


0.17 


2.87 


Power law plus 


















absorption 


1.51/23 


2.27 ±0.19 


2.78 ±2.97 


12.98 


6.49 


5.96 


0.04 


2.75 


F = 2 power law 


















plus absorption 


1.57/24 


2.0 


0.04 ±1.66 


15.0 


6.19 


7.92 


0.14 


2.68 


Power law plus 


















fixed absorption 


1.48/24 


2.18±0.11 


1.0 


13.38 


6.36 


6.40 


0.10 


2.82 



All fits include a Galactic column of 3.2x10^° cm 2. Quantities in bold are fixed during fitting. Fluxes are 
quoted in the observed frame. Luminosities are quoted in the rest frame. 
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Fig. 1.— EPIC pn spectrum of SDSS J1030+0524, 
plus a single power law fit & fit residuals with 
absorption fixed at the Galactic value. The MOS 
spectra have not been plotted for clarity. 




Hydrogen Column 

Fig. 2. — Confidence plot of F vs. -/Vh (units are 
10^2 cm~2) from the fit where both parameters are 
allowed to vary. Contours are la, 2a & 3a. 



